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Abstract 

A transfer-matrix algorithm is presented herein as a beginning to study the transmission charac- 
teristics of coherent hght through three-dimensional periodic microstructures, in which the struc- 
tures are treated as two-dimensional-layer stacks and multiple reflections are considered negligible. 
The spatial-correlated noise is further introduced layer by layer to realize the actual decoherence 
of the light and allows for statistical investigation of the partial spatially coherent optics in trans- 
parent mediums. Numerical analyses show comparable results to the Gaussian Schell model in 
free-space cases, indicating the validity of the algorithms. 

PACS numbers: 42.25.Kb, 78.35. +c, 42.25.Fx, 05.10.Gg 



* Electronic address: lsunta.ho@msa.hinet.netl 



1 



I. INTRODUCTION 



Nanoscale structures have achieved novel functions in electro-optic devices such as optical 
filters, optical modulators, phase conjugated systems, optical attenuators, beam amplifiers, 
tunable lasers, holographic data storage and even as parts for optical logic systems jl- 14| 
over the last few decades. In contrast to the assumption of the ideal coherent (plane-wave) 



light in most theoretical studies 15l-l2l|. the associated decoherence characteristics and in- 
tensity distribution of the light, e.g. by light-emitting diodes (LED) 22l-24|, otherwise 
introduce another critical issue for both fundamental research and potential industrial ap- 
plications. Motivated by these concerns and based on our previous study on coupled-wave 



theory 



20l |. we present a transfer-matrix algorithm as a beginning to study the propagation 



of coherent light through three-dimensional periodic microstructures. This work discretizes 
the structures into two-dimensional-layer stacks composed of isotropic or birefringent mate- 
rials in an arbitrary order, and multiple reflections are considered negligible. By employing 
a concept similar to the Langevin dynamics in the finite difference method 25|, a noise 
phase characterized with spatial correlations r is further included layer by layer to mimic 
the actual decoherence of the light and allows for statistical investigation of the partial spa- 
tially coherent optics. Here, the correlation length r of the noise is assumed to be larger 
than the wavelength of the incidence A, because the intense decoherence (small r) causes 
strong scattering and the multiple reflections thereby cannot be ignored. The generation of 
the noise function is fulfilled by using a standard numerical technique based on the Fourier 



transform [26|, l27| and is formulated in the appendix for the two-dimensional cases studied. 
Numerical analyses for the propagation of the free-space Gaussian beam and the diffractions 
of the beam incident through a liquid- crystal gra ting at different spatial decoherence show 
comparable results to the Gaussian Schell model 



28 



29| and can validate this study. 



II. THEORETICAL FORMULAE 



This section neglects the multiple reflections and derives the transfer-matrix algorithm 
that is much easier to manipulate algebraically, yet accounts for the effects of the Fresnel 
refraction and the single reflection at the interfaces of the medium. The assumption of no 
multiple reflection is legitimate for most practical transparent materials. In what follows, a 
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FIG. 1: Geometry of three-dimensional RCWA algorithm for a multi-layer stack with two- 
dimensional periodic microstructures in arbitrary arrangement with isotropic and birefringent ma- 
terials . 

reformulation, including the noise phase of electromagnetic fields, is fulfilled to demonstrate 
the decoherence behaviors of the light. 

A. Transfer matrix method by RCWA 

Referring to Figure [H the structure is stratified along the medium normal into layers 
and these layers contribute to a total thickness of the stacks = Xl^Li ^i- Each layer is 
considered as an arbitrary arrangement with homogeneous or birefringent materials, and the 
dielectric coefficients are generally described by a matrix: 
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with 

£xx = nl+ (nl - nl) sin^ 9o cos^ 0o, 
£xy = £yx = {nl - nl) sin^ 60 sin (po cos 
£0:2 = £zx = {nl - nl) sin 9o cos cos 0o, 

= n^+(n^-n^)sin^^„sin^0„ 
£2/2 = ^zy = ('T'e ~ ^o) sin COS sin (f)o, 

Szz = (n^-n^)cos2^<„ (2) 

Here, Ug and are extraordinary and ordinary indices of refraction of the uniaxially bire- 
fringent medium, respectively, 9o is the angle between the optic axes and the z axis, and (po 
is the angle between the projection of the optic axes on the xy plane and x axis. 

To introduce the rigorous coupled-wave theory to the stack, all of the layers are defined 
as having the same periodicity: A^- along the x direction and Ay along the y direction. The 
periodic permittivity of an individual layer £ in the stack thereby can be expanded in the 
Fourier series of the spatial harmonics as: 

£ij (x, y, zi) = ^ Ei^^gh (ze) exp + (3) 

g,h ^ ^ 

Here, we have defined variables ko — ui^iiqEq = X' ^ ~ ^ ~ \f%'' ^ ~ ^ ~ 
y — koy, and z — koz. X is the vacuum wavelength of the incident wave. Note that £ije{x,y,z} 
arc defined as functions of position {x,y,z) and Eij are defined for {x,y,z). A parallel 
transform for the electromagnetic fields through the stack is expressed in terms of Rayleigh 
expansions: 

V^^ix,y,ze) = '^egh{ze)exjp[-i{nxgx + nyhy)] (5) 

9,h 

\/%H(x,y,z^) = '^hgh{ze)eyip[-i{nxgx + nyhy)] (6) 
f^xg = ni sinO COS (j) — g— (7) 

A-x 

nyh = nism.9s\D.(j) — h— (8) 

Ky 

This indicates the propagation of light along the direction ngh = n^gi + ny^j + $,ghk with 
^gh — {njf^E) ~nyhnyh—nxgnxg)^^'^ , in which nj and ue are the refraction index and correspond 



to the propagation in the incident and emitted regions, respectively. 9, are the incident 
angles defined by sphere coordinates, and z is the normal direction for the xy plane of 
periodic structures. 



Organizing the algorithms from our previous study 



20| while ignoring the multiple re- 



fiections, the transfer-matrix formulae for the microstructure can thus be written as: 
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(9) 



Here, Mjgji^Tv} is the transfer matrix corresponding to the ith (anisotropic) structured layer. 
It is formulated by the eigen- vector matrix T with column eigen- vectors and the diagonal 
eigen- value matrix /t of the characteristic matrix by: 



^x^zz^zy 
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(10) 



(11) 



where the notation (.) represents Ng x vectors with components (.)gh, rix ijiy) are NgN^ x 
NgNh diagonal matrices with the diagonal elements n^g (nyh), and eij(z^x,y,z} are NgNhxNgNh 



matrices with elements eij„h formulated by / ~ eijf or /, 



Igh 



Yliq'h' ' 



'ij,(g-9')ih-h')Jg' 



h')fq'h'- 



^g{h) define the number of considered Fourier orders g (h) in the x (y) direction. 1 is 
NgNh X NgNh identity matrix. One may understood that (T.'^'')"^ term represents the 



coordinate transformation from the components of tangential fields 



^x,i hy^i dy^i hx^i 



at ith interface in Equations (EJ-dS]) into the independent components of the eigen- modes of 
the ith layer, exp [iHi^i] term then follows the corresponding eigen-mode transition lasting a 
distance ^j, and T^^"'' term is the inversely coordinate transformation back to the components 
of the tangential fields defining ifi^i at the next interface. The superscript t indicates the 
matrix transposition. Summarily, Mj demonstrates the propagation of fields (f^j to iii-^i) 
through the ith layer. 

In the (isotropic) uniform incident {i = 0) and emitted {i = N + 1) regions, the eigen- 
modes are chosen as Eg and Mg {E' and Mq) |20|, l2l|, thereby representing the physical 



forward (backward) TE and TM waves, i.e. transverse electric and transverse magnetic fields 
corresponding to the plane of the diffraction wave, respectively. The transform between the 
eigen- modes and the tangential fields f^ Q = [cxfi hy^ Cy^ /i^.o]* in the incident region {i = 0) 
is written as: 
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(12) 



Here, n„ and are NgNh x NgNh diagonal matrices with normalized elements and 
respectively. is the matrix with elements l/^gh (not the inverse of the matrix ^), in which 
mgh = {nyhUyh + nxgnxgY/"^, ^gh = {s I " UyhUyh - rixgnxgY^^, and ej = n] have been defined 
for the incident region. A similar transform for j^^^ in the emitted region can be derived 
straightforwardly by replacing all the Si in Equation ( 1T2|) with se and can be obtained as 
%N+i = T^el[E+^+i M+jv+i ^9>+i Mg,N+i\^ in which = {ee - riyhUyh - rixgnxgY^^, and 
Ee = n\. 

Moreover, Me„t is the matrix representing the light propagation from the incident region 
into the medium, and indicates the essential refraction and the reflection at the first interface 
of the medium. To consider these effects in a simple way, a virtual (isotropic) uniform 
layer, which has zero thickness and (scalar) average dielectric coefficient Ea = n^^g, e.g. 
navg = (^e + ^o)/2 for the liquid- crystal grating, is assumed to exist between the incident 
region and the Igt layer. Me^t thereby can be obtained as: 

Wi"^ 
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ent 



Wi W2 
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(13) 
(14) 



Here, Ti*J is defined in Equation ( fT2l) with the replacements of Ej by Ea, 6 



gh 



^yh^yh ^xg^: 
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Parallel to the argument of Ment, another similar 
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virtual (isotropic) uniform layer exists between the emitted region and the Nst layer, and 
M.ext is approximated as: 



M 



ext 







(15) 
(16) 



B. Algorithms extended for optical decoherences 

To demonstrate the spatial decoherence of the light, the electric fields of an incident 
unit-amplitude plane wave are first denoted as: 



E (r) = Eo exp [—ik ■ r + ir] (r)] 



(17) 



The time dependence of exp (iut) is assumed and omitted here. ri{x,y, z) is the noise phase 
of E and is conditioned by the correlation relation with a characteristic spatial coherent 
length r: 



(r/(r)r/(r')) = 7 (r - r) 



TT exp 



2r2 



(18) 



with {r] (r)) = 0. The correlation function 7 (r — r') is chosen as Gaussian distribution and 
describes a strong (weak) phase correlation between fields of closer (farther) positions than 
r. In the limit of r — )■ 00, the 77 (r) is constant and the E field in Equation ffTTl) shows 
the coherent incident wave, whereas in the limit of r — )■ 0, an incoherent (total random 
phase) wave is exhibited. Here, the generation of the noise function 77 (r) is fulfilled by using 
a standard numerical technique based on the Fourier transform and is formulated in the 
appendix for the two-dimensional cases studied. To realize the spatial decoherence analyses 
for the light, we introduce the concept similar to the Langevin dynamics for the stratified 
layers along the z direction and include the noise phase layer by layer in Equation ([9]): 
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(19) 



Here, R indicates the operations of three numerical processes on ii,ie{ir^N+i} = 
Mi_i...RM2RMiRMe„f ■ [^+o M^,o E^fl ^^gTo]* at interface by: (a) first fulfiUing the 
inverse Fourier transform from the k components of the tangential fields f£j(x,|/) into the r 
components, i.e. {^/YqEi, a/ZqH ) by Equation (111); (b) next including the noise phase by 
multiplying the spatial r components of the fields by the noise phase e*'''-^'^'^-*, or e*''*^^'^'^-' 
in our work, as Equation ( |T71) : (c) and finally performing the Fourier transform on the 
r components to obtain the k components of the tangential fields if:^{x,y) including the 
noise phase. Hence, the field jg{i^jv+i} through the medium can approximately mimic 
decoherence behaviors by these processes. 



III. NUMERICAL RESULTS FOR TWO-DIMENSIONAL EXAMPLES 

To study the decoherence behavior of light, we first analyze the propagation of the Gaus- 
sian beam with different coherent lengths in a two-dimensional xz free space, and fulfill a 
statistical numerical analysis over an ensemble of systems with 500 identical iterations. The 
characteristic period of the structure (A^: = 48/im) is chosen to be much larger than the 
Gaussian beam profile (~ e-^^/2o-^ with a = 6fim) so that the periodic beam can be consid- 
ered to be isolated in near-field optic analyses. The relevant parameters are given as: the 
x-grid size A^, = 0.05/im, the thickness of the stacked layer Zi^^i^^} = 0.05/im, the wave- 
length A = 0.55/im, the y-polarized incidence y/YoEy{x, z = 0) = e~^'^/^'^^, the refraction 
indices ni = ue = rio = 1.0 ~ = 1.0 -|- 10~^, and the Fourier order —50 < (7 < 50 (/i = 0). 
Figure [2] illustrates the intensity patterns \^/YQEy{x^ + \ y/ZQHx{x, z)p of the free-space 
Gaussian beam with r = 5;um, r = l/im, and r = 0.5/im along the light propagation. Figure 
|3] shows the corresponding numerical intensity to the Figure [2] after propagating a distance 
at 2; = 35yum. The results indicate that the stronger decoherence (smaller r) of the light 
causes more a intense scattering along the propagation and shows comparable conclusion 
from the Gaussian Schell model 



281. 

To further investigate the decoherence behaviors of the light through the microstructures, 
a liquid-crystal grating film is included to introduce the diffraction of light. Figure H] shows 
the one-period orientations of the liquid-crystals directors in a single-layer film. For this 
case, rather than employing the definition of the grating period h^x^ic = 1.2/xm as in Figure 
m we apply SRC-RCWA scheme and set the characteristic period A^; = 48/im = 24 x K^^lc- 
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(a) (b) (c) 

FIG. 2: Intensity patterns of the free-space Gaussian beam with coherent length (a) r = 5^m, (b) 
r = 1/im, and (c) r = 0.5/xm for the ensemble of 500 iterations in the 48/im x 40/Ltm xz plane. 

1=5.0 flm 




10 20 30 40 
x(nm) 



FIG. 3: Numerical intensity of the Gaussian beam corresponding to the results of Figure [2] with 
r = 5^m, T = l^m, and r = O.bfim after propagating a distance at z = 35fim. 

Hence, the long-range profile of the Gaussian beam in spatial space can be considered, and 
simultaneously the small-angle emittance of the beam (~ l/A^) defined by the Fourier com- 
ponent in Equations (C])-® can be demonstrated. Figure [5] illustrates diffraction patterns 
\\^Ey{x,z)\'^ + \ y^Z^Hx{x, z)]"^ of the beam incident through a forty-period hquid-crystal 
grating (enclosed by dashed lines) in 48/im x AOfim xz plane, in which the coherent lengths 
of the beam are statistically considered to be (a) r = 5yum, (b) r = 1/im, and (c) r = 0.5/im 
for an ensemble of 500 iterations. Figure [6] shows the corresponding numerical intensity of 
the beam to Figure |5] after propagating a distance aX z = 35fim. The results indicate that 
the coherent incidence (r = 5/im) through gratings leads to definite diffractions as described 
in most studies |20i], while the decoherent one (r = O.S/xm) exhibits an intense scattering 
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FIG. 5: Intensity patterns of the Gaussian beam incident through a forty-period grating (dashed- 
hne rectangle) in the 4:8fxm x 40/im xz plane. The coherent lengths of the beam are statistically 
considered to be (a) r = 5/im, (b) r = Ifim, and (c) r = 0.5^m in an ensemble of 500 iterations. 

similar to that in Figure M^c) and Figure [3l The incidence with r = l/im in Figure Et^b) 
and Figure [6] shows a simultaneous diffraction and scattering results and deviates from the 
results by the coherent optics. 

IV. CONCLUSIONS 

This work has presented the formulas of the transfer-matrix method to conduct deco- 
herence analyses in three-dimensional periodic microstructures. The algorithms are also 
devoted to doing a per-study of the light through turbid mediums that are related to the 
interaction between the fluctuated electron/dipole motions and the decoherence behaviors 
of light. Two numerical analyses for the propagation of free-space Gaussian beams and the 
diffraction of the liquid-crystal grating are then applied to verify the validity of this work, 
obtaining reasonable results. 
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FIG. 6: Numerical intensity of the Gaussian beam corresponding to the results of Figure [5] with 
T = 5^m, T = Ifim, and r = 0.5fim after propagating a distance at z = 35fim. 
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Appendix A: Generation of two-dimensional spatial correlated noises 

In the following context, we introduce the generation of fluctuation function r] (r) by 
the discrete Fourier method, and thereby realizing the E field describing partial coherence 
(0 < r < oo) in Equation (fTTll . First, we introduce //(k) in 2D momentum space, which 
follows the transforms of: 



or alternatively the discrete representation in N x N grid data under spatial interval A for 
our case 




(Al) 



(A2) 
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N-1 



Vg,h 



^ exp 



n,t)=0 

N-1 



'ij;^{g-u + h-v) 



g,h = Q,l,...,N -I (A3) 



= 0,1,...,A^- 1 (A4) 



Here m, f (fifj/i) are the indices for the spatial (momentum) space. The correlation function 
of the momentum space corresponding to that of the spatial space in Equation (fT8|) is next 
derived as: 



(r7(k)r7(k')) = 47r'5 (k + k') 7 (k) 

= 47r^(5 (k + k') ■ 27rV^ exp 



1 2 2 

-k r 



(A5) 



where 7 (k) is the Fourier transform of 7 (r). The discrete representation of Equation (1A5I1 



IS 



{Vg,hr]g',h') = N'^^^lg,Jg,N-g'5h,N-h' (A6) 

Here, 7g,/t is the discrete representation for 7 (k) and exhibits the symmetry of 7g^/i = 
lN-g,h = 1g,N-h = lN-g,N~-h- Accordingly, by the discrete correlation equation in Equa- 
tion flA6p . the fluctuation function rjg^h can be generated in the discrete momentum space 
by: 



r]g^h = NA,/y^ag^h (A7) 

in which ag^h is the Gaussian random number (complex number) with the average being 
zero, and it is conditioned as (ag^hC^g'y) = Sg^iy^g'6h,N-h', which corresponds to Equation 
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(lA6p . Definitely, these requirements for ^ can be realized by a simple process as: 
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<^N-g,N-h = 





0-g,N-h + ibg^N-h) , OtN-g,h — —f^ {o-g,N-h — ibg^N-h) (A8) 

Here, a^^/j and 6g ^ are independent Gaussian random numbers (real numbers) with an av- 
erage of zero and variance of one. Note that it is necessary to generate N x N independent 
Gaussian random numbers for N x N fluctuations of r^^ /j. 

Finally, ?7„^„ in the discrete spatial space can be obtained straightforwardly by the 
Fourier transform of rjg^h in Equation IA41 such that the electric field E (x, z) = 
Eoexp [—ik ■ r + r] {x, z)] in the studied xz region can be evaluated. A further treatment 
of rju^v Vu,v — V, in which fj means the average noise over the N x N grid data, is executed 
to ensure the condition {f]{r)) = 0. It is straightforward to extend to arbitrary N x M 
periodic grid data. Numerical results of the noise phase with r = 5,1, 0.5/im are illustrated 
in Figure [7] for reference. 
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